Propagation of a Dark Soliton in a Disordered Bose-Einstein Condensate 
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We consider the propagation of a dark soliton in a quasi ID Bose-Einstein condensate in presence 
of a random potential. This configuration involves nonlinear effects and disorder, and we argue that, 
contrarily to the study of stationary transmission coefficients through a nonlinear disordered slab, 
it is a well defined problem. It is found that a dark soliton decays algebraically, over a characteristic 
length which is independent of its initial velocity, and much larger than both the healing length and 
the ID scattering length of the system. We also determine the characteristic decay time. 
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Phase coherent systems display wave mechanical pro- 
perties distinct from those typically observed at macro- 
scopic scale. In particular, transport in presence of dis- 
order is strongly affected by interference effects, leading 
to weak or strong localization, as observed in many dif- 
ferent fields (electronic or atomic physics, acoustics, elec- 
tromagnetism). Our understanding of these effects have 
made great progresses over the last decades in the case of 
non-interacting linear waves. Some studies have conside- 
red the propagation of plane waves or bright solitons in 
a disordered region in the case of attractive interaction 
(see, e.g., the review pj and the discussion below of the 
results of Ref. 0), but almost nothing is known in the 
case of repulsive nonlinearity. 

The field of Bosc condensed atomic vapors allows new 
investigations of such phenomena in presence of repul- 
sive or attractive interaction, in an intrinsically phase 
coherent system over which the experimental control is 
rapidly progressing. Such studies have begun with the 
observation of "fragmentation of the condensate" over a 
microchip random potentials have recently been en- 
gineered using an optical speckle pattern p|; and it has 
also been proposed to implement disorder by using two 
different atomic species in an optical lattice |5|. 

In the present Letter, we study the transport proper- 
ties of a quasi one dimensional (ID) Bose-Einstein con- 
densate in presence of disorder and repulsive interaction. 
The configuration we study corresponds to the "ID mean 
field regime" 0, where the system is described by a ID 
order parameter ip(x,t) depending on a single spatial 
variable: the coordinate x along the direction of propa- 
gation. ip(x,t) obeys the non linear Schrodinger equation 

where U(x) is the random potential and g an effective 
coupling constant which reads g = 2Hu)±a in the case 
of particles experiencing an effective repulsive interac- 
tion characterized by the 3D s-wave scattering length a 
(a > 0), and a transverse harmonic confinement with 
pulsation w± Q. It is customary to define the oscilla- 



tor length a±_ — (ft/mwi) 1 / 2 and a% = a 2 L /2a (— a\ is 
the ID scattering length Denoting by n%o a typical 
value of \4>(x,t)\ 2 , the ID mean field regime corresponds 
to a situation where 1 <C n\ D a\ <C (ai/a±) 2 . The first 
inequality ensures that the system does not get in the 
Tonks-Girardeau limit and the second that the transverse 
wave- function is the ground state of the linear transverse 
Hamiltonian [alS]- 

A particular issue specific to Eq. (JTJ is the very pos- 
sibility to define a transmission coefficient. Since the 
equation is nonlinear, it is not possible in general to dis- 
entangle an incident and a reflected current in the re- 
gion upstream the potential (in other words, a reflected 
atom will interact with the incident beam) . A possible 
way for avoiding this problem is to change the transverse 
confining potential upstream the disordered potential so 
that, in this region, nonlinear effects become negligible 
|9|. However, even in this technical difficulty 

arises because of multistability: several stationary solu- 
tions exist for a given input state Q, Q. Moreover, in 
the case of repulsive interaction we consider here, for ex- 
tended enough disordered region, no stationary solution 
exist and the transmission coefficient can only be defined 
via a time average [UJ . 

A way out of these difficulties consists in studying the 
propagation of a soliton in the system. This constitutes 
an intrinsically time-dependent problem, but the input 
and output states can be precisely characterized, and the 
transmission is simply defined by comparing the large 
time behaviors (t — > ±oo) of the solution. This route 
has been followed by Kivshar et al. Q in the case of 
attractive nonlinearity (g < 0). In this solitary 
wave is a bright soliton, characterized by two parameters: 
its velocity V and the number of particles N inside the 
soliton. The disordered potential was taken as 

h 2 

U (x) = g imp Y, s ( x - x n) , where g imp = — . (2) 
* — ' mo 

n 

U(x) describes a series of static impurities with equal 
intensity and random positions x n . The x n 's are mi- 
correlated and uniformly distributed with mean density 
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n imp . In this case <U(x)>= g imp n imp and <U(x{)U(x2)> 
- <U(x\)> x <U(x 2 )>= (h 2 /m) 2 D 8(xi - x 2 ), with 
D = n imp /6 2 . From what is known in the case of linear 
waves, this type of potential is typical insofar as localiza- 
tion properties are concerned. In the weakly nonlinear 
regime mV 2 /2 > h 2 N 2 /(mal), it was found in Ref. Q 
that the soliton velocity remains approximatively con- 
stant in the disordered region, whereas N shows an ex- 
ponential decay similar to what occurs for a linear wave 
packet. In the opposite strongly nonlinear regime, it was 
found that the soliton behaves very differently, leading 
asymptotically to a configuration where both N and V 
become practically constant (independent of x). 

In the present Letter, we consider the case of repulsive 
nonlinearity where the solitary waves are dark solitons. 
We find that the propagation of these solitons in a dis- 
ordered potential is quite peculiar for the two following 
reasons: first, the strongly and weakly nonlinear cases 
cannot be considered as distinct because, in a given sys- 
tem, the number of particles forming the soliton cannot 
evolve independently from its velocity; and secondly, a 
dark soliton has a velocity bounded by the velocity of 
sound in the system. As a result, dark solitons behave 
differently from the bright ones studied in Ref. 0]: ini- 
tially rather "nonlinear" solitons decay algebraically (and 
not exponentially) , becoming eventually "linear" . Be- 
sides the length covered by the soliton in the disordered 
region is independent of its initial velocity. 
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c = h/m!;. A dark soliton has a velocity V < c, an energy 
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Figure 1 : Density profile of a dark soliton incident with 
velocity V = c/2 on point-like repulsive obstacles with 
random positions corresponding to a potential U(x) given 
in Eq. J2J), with £/b = 0.2 and n imp £ ~ 0.1. TJie arrow 
represents the direction of propagation of the soliton. 

Let us thus consider a dark soliton with initial velocity 
V, incident from the left on a disordered potential of type 
with xo = < x\ < .... This situation is illustrated 
in Fig. 1. The soliton is characterized by two parame- 
ters, its velocity V and the asymptotic background den- 
sity rioo = lim^^ioo \ip(x,t)\ . Instead of rioo, one can 
equivalently employ the chemical potential fi = g Hoc , 
the healing length £ = (ai/rioc) 1 / 2 , or the speed of sound 
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and consists in a density trough of typical extension 
£ (1 — y 2 /c 2 ) -1 / 2 , corresponding to a number of missing 
particles AN = 2 (ai/£)(l - V 2 /c 2 ) 1 / 2 . In the ID mean 
field regime where Q holds, ai ^> £ and AN is typically 
a large number, except in the limit where V is close to c. 
This occurs at velocities around V ctit = c [1 — (^/2ai) 2 ] 1 / 2 . 
At such velocities AA^ ~ 1 and the soliton has an exten- 
sion ~ Ol. 

We consider the case where the average separation 
between the impurities is much larger than the healing 
length (n imp £ <§; 1) and the initial velocity of the soliton 
is not close to c. In this case, the scattering of the soliton 
from the impurities can be treated as a sequence of inde- 
pendent events. When the soliton encounters a single ob- 
stacle, it radiates phonons which form two counter prop- 
agating wave packets moving at velocity c. Accordingly, 
its energy decreases by an amount 8E sol — —E+ d — E~ d , 
where E+ d (E~ d ) is the forward (backward) emitted en- 
ergy. It was found in Ref. 0] that 



Mr) F± (V/c), 



(4) 



where i ?± (u) is a dimensionless function defined for v — 
V/c S [0, 1] as 

,+oo y 4 (-v± v/l + y 2 /4) 2 
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Equation Q is a perturbative result valid in the limit 
b ^> £ and V 2 3> c 2 (£/b). The first inequality ensures 
that the impurity only weakly perturbs the static back- 
ground and the second that the scattering of the soliton 
by the impurity can be treated perturbatively. The soli- 
ton having lost energy during the collision, its velocity 
changes by an amount SV = cSv which, from |J5J, is re- 
lated to 5E eol via vSv (1 - v 2 )- 1 = -\SE e J E eoX . 

Since n imp £ <C 1, one can go to the continuous limit 
considering the successive collisions as a sequence of ran- 
dom uncorrelated events. Over a length 6x the solitons 
will experience n imp i5a; such collisions. This leads to the 
following differential equation: 



dv _ 1 F+(v) + F~(v) 
dx 4 .t v \/l - v 2 



(6) 



where xq = ai& 2 /(£ 3 n imp ) = ai/(Z?£; 3 ). Equation © can 
be solved analytically in the high velocity regime, when 
v -> 1. In this limit, F + (v)+F~{v) ~ 4(l-« 2 ) 5/2 and, 
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for a soliton of initial velocity V i: 

V(x) f 1 l-(Wc) 2 



= <1- 



one obtains 
, 1/2 



l + [l-(Wc) 2 



2x 
15 i'n 



(7) 



We compare in Figure 2 the results of this approximate 
solution with the numerical solution of Eq. © in the 
cases Knit/c = 0.75, 0.5 and 0.25. The agreement is very 
good, even for initial velocities which are not close to c. 
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Figure 2: Evolution of the velocity V of a dark soliton 
as a function of the distance x traveled in the disordered 
region. In each plot, the thick line corresponds to the 
numerical solution ofEq. JjJ) and the thin line to Eq. Q). 
Case (a) corresponds to an initial velocity V m — 3 c/ '4, 
case (b) to V Mt = c/2 and case (c) to V m — c/4 . 

The soliton is accelerated as it progresses through the 
disordered region (as seen in Figure 2) because it radiates 
energy at each collision with an impurity. This increased 
velocity after a loss of energy is a typical feature of dark 
solitons which can be considered as particles with a nega- 
tive kinetic mass which decreases with increasing energy 
[l2T| (see Eq. ©). One also notices in Figure 2 that V 
saturates when it gets close to c, meaning that, in this 
regime, the rate of energy loss decreases. The reason 
for this phenomenon is that a dark soliton cannot have 
a velocity higher than c. As a result, when its veloc- 
ity reaches this upper bound, the soliton cannot lose a 
large fraction of its energy, because this would lead, after 
the collision, to an unphysical value of V (larger than c). 
This phenomenon has an important consequence on the 
maximum distance L over which the soliton can travel 
in the disordered region. As seen on Figure 2, L is very 
large and seems independent from the initial velocity of 
the soliton. In order to get a quantitative evaluation, we 
define L as being the length after which the soliton is 



a trough containing only one particle, i.e., the velocity 
V(L) in Eq. (Q) reaches the value V ctit . In this limit the 
soliton can no longer be detected by standard imaging 
techniques, and for all practical purposes one can con- 
sider that it has totally decayed. From J7J) one obtains 
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This confirms what was inferred from Figure 2: a slow 
soliton will initially decay more rapidly than a fast one 
and altogether, the distance over which solitons can 
travel before completely decaying is independent of their 
initial velocity. As expected, L decreases for increasing 
disorder, the effect of the disorder being measured by the 
dimensionless parameter D £ 3 , i.e., by the two points cor- 
relation function of the random potential. Irrespectively 
of the value of the parameter D£ 3 , we remark that L 
is large compared with ax, since in the ID mean field 
regime a\ £. Hence, a dark soliton covers quite a large 
distance in the disordered region before decaying. 

The distance L is covered in a time r which we now 
evaluate. It is important to realize that V is not the 
average velocity of the soliton, but its velocity between 
two obstacles: in vicinity of an impurity, the velocity of 
the soliton decreases if g imp > and increases if g imp < 
0. As a result, the asymptotic position of the soliton 
is shifted compared to what it would be in absence of 
obstacle. In the case of a single impurity, this shift A can 
be quite accurately evaluated by means of the "effective 
potential theory" as being [ll"| 



A = 



d.v 



1 - 



1 



yjl - U c!f (x)/mV 2 



(9) 



where {/„« is an effective potential which reads in the 
case of a point like impurity U o{f (x) = cosh -2 (x/£). 
In the limit V 2 > c 2 (£/6) where Eq. |gj holds, the 
shift reads A ~ — c 2 ^ 2 /(2bV 2 ). In the presence of mul- 
tiple impurities, going to the continuous limit, one ob- 
tains that during a time St the soliton covers a distance 
5x = VSt + (n imp 5x) A. Combining this relation with 
Eq. © one obtains a differential equation allowing to 
determine v — V/c as a function of t: 

dv _ c F+(v) + F~(v) 
dt Axq 



(10) 



In the limit v 
solution 
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1, this equation admits the analytical 
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We compared this approximate result with the numerical 
solution of iJTUj l where A was evaluated through and 
found that the accuracy of ltTT|) is always very good, even 
for initial velocities not close to c (as was also the case 
for the approximate expression Q). 

The decay time r of the soliton is the time at which 
V = y„ it / C ~l-i(£/2 ai ) 2 : 

r = -(l + ^f). (13) 

In this expression - as in JSJ - we neglected a corrective 
term depending of the initial velocity, smaller by a fac- 
tor (£/ai) 2 than the leading term, r is proportional to 
L/c, with a slight modification due to the shift induced 
at each scattering 0|: repulsive obstacles (b > 0) lead to 
an increased decay time since the soliton covers the dis- 
tance L slightly more slowly that in the case of attractive 
obstacles. 

Formula (|13[) can be given a simple physical interpre- 
tation (in a less rigorous setting) in the framework of 
the "effective potential approximation" ^lj. I n this ap- 
proximation, solitons are considered as classical particles 
of mass 2m evolving in a potential U etl . One thus has 
<mx 2 + U c!f (x) >=< mV 2 (x) >. The mean value of 
U ctt is the same as the one of U 0] and from Figure 
2, one sees that at leading order it is sensible to ap- 
proximate <mV 2 (x)> by mc 2 = /i. One thus obtains 
<i 2 >~ c 2 (l— <U(x)> Iii). Finally, r can be evaluated 
through the formula 

L L Lf <U(x)> \ 

T ~ <x> ~ <i 2 >V2 - c + 2f i J 1 1 ' 

which is identical to . Since formulas © and i|14|) de- 
pend only on simple characteristics of the random poten- 
tial (the average and the two points correlation function) , 
we expect them to be of very general validity, poorly af- 
fected by the specific potential present in the disordered 
region. 

A final point to clarify is the effect of the random 
potential on the occurrence of superfluidity and Bosc- 
Einstein condensation, i.e., is Eq. truly applicable ? 
In the strong disorder limit, a quantum phase transition 
occurs at T = leading to a (non superfluid) Bose glass 
phase 15] where the description of the system with a sin- 
gle order parameter ip(x, t) is inappropriate. However, in 
the case we consider here of an atomic vapor described 
as a weakly interacting Bose gas, it has been shown that 
a small amount of disorder does not drastically alter the 
properties of the system, but merel y d ecreases the con- 
densate and the superfluid fraction [lj|. More precisely, 
based on the evaluations presented in Ref. 01 one can 
show that this effect is negligible provided n ilnp ^ 3 /b 2 <C 1, 
which is the case in the present study. 

In conclusion, we have presented a description of the 
motion of a dark soliton in a disordered region. The 



soliton radiates energy when it encounters an obstacle. 
The repulsion between the particles has important con- 
sequences on the propagation of the dark soliton, whose 
salient features are all at variance with the one expected 
in the case of a linear wave packet or of a bright soliton: 
(i) the soliton is accelerated to the velocity of sound and 
disappears, (ii) its decay is algebraic, and (iii) the char- 
acteristic decay length and decay time are independent 
of the initial velocity of the soliton. 

These results are generic and apply to many different 
fields (among which, optics in nonlinear fibers with pos- 
itive group velocity dispersion) but the most promissing 
experimental relizations seem to be achievable for a Bose 
condensed atomic vapor, either in a corrugated magnetic 
guide over a microchip 3], or in an elongated trap in 
presence of an optical speckle pattern |4| . 
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